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Abstract 

A new method to efficiently downdate an estimate and covariance generated by a discrete time 
Square Root Information Filter (SRIF) is presented. The method combines the QR factor down- 
dating algorithm of Gill [3] and the decentralized SRIF algorithm of Bierman [4]. Efficient removal 
of either measurements or a priori information is possible without loss of numerical integrity. More- 
over, the method includes features for detecting potential numerical degradation. Performance on 
a 300 parameter system with 5800 data points shows that the method can be used in real time and 
hence is a promising tool for interactive data analysis. Additionally, updating a time- varying SRIF 
filter with either additional measurements or a priori information proceeds analogously. 


Introduction 

A typical 24 hour data arc for the GPS demonstration on TOPEX/POSEIDON [6] will contain 
« 30,000 data points. To process these data points on a VAX 8530, a sequential SRIF filter will 
require nearly 3 hours of CPU time. Upon processing, should outliers be discovered in the data, due 
to cycle slips, atmospheric fluctuations, multipath, etc . . ., they must be removed from the data arc 
and another 3 hours of CPU time would then be needed for reprocessing. The method presented 
here permits efficient removal of these outlying data points without reprocessing the entire data 
set. 


The discussion is organized in the following manner. First, a decentralized approach to up- 
dating a time- varying SRIF with a single measurement is presented. From this development, it 
is evident how to remove a measurement from a time-varying SRIF. Next, the one-component- 
at-a-time process noise methodology is applied in the actual implementation. The necessary linear 
combination of the data equation and the smoothing coefficients is then presented. Finally, the 
efficiency and the numerical integrity of the method is discussed. 

Updating a Tune- Varying SRIF with a Single Measurement 

The goal is to find estimates (aro, , • . . , to minimize the the least-squares performance 
functional 
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for the linear discrete-time dynamic model 


with measurements 


and the single measurement 
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where () (j) represents a priori information between time updates j and j+1. The data noise Vj 
and process noise w 3 are independent mean zero white Gaussian noise processes with covariances 
I and R w (j) -1 respectively. Of note, when a pseudo-epoch state formulation is used <f)j 

and G have a simplified structure. 
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Ignoring the additional data equation yk = Hk%k + Vk> the SRIF alternates between measure- 
ment updates 
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f °7 ~ °’ Ar “ L The e * term ^Presents the residual sum of squares from the i’th mesurement 

update; the ()* terms represent the smoothing coefficients. Additionally, the notation f) . represents 
quantities that include data up to time j 1. ^ 

To update the filter with the additional data equation y k = H k x k + i/ t) it is necessary to 
express this equation in terms of x A+1 and w k . 
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This equation is then merged with the smothing coefficients of the SRIF. 
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The ()** terms replace the old smoothing coefficients and y *+ 1 = H k + 1 ** +1 + is the ensuing 

additional data equation. The time updates in equation (10) and the merges in equations (11) 
and (12) are then repeated for k + 2, k + 3, . . . , N. In the end, the terminal data equation y N = 
Hn*n + UN is then merged with the terminal SRJF array. 
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The notation () w represents quantities resulting from merging the terminal data equation with the 
terminal SRIF array. The filter estimate and covariance are then R^ z N and R N R N .respectively. 
To obtain the smoothed estimates and covariances, the Dyer-McReynolds smoother [2] would then 
operate with the ()* terms for j = 0, k - 1 and with the ()** terms for j = k,N- 1. 

Downdating a Time- Varying SRIF with a Single Measurement 

The last lines of equations (11) and (12) contain the key to reversing this process If the 
data equation y k = H k x k + u k has been included in the SRIF, then the ()** terms of (12) are 
available. Additionally, the data equation y*(k) = H*,{k)w k + H*, x (k + l)ar x +i + ”k in l 11 ) can 
be generated from the measurement that is to be removed. It is possible, as shown below, to solve 
for a transformation to produce the upper ()* terms of (11) and the ensuing data equation y*+i - 
R* + i** + i +u k+ 1 in (12). The upper ()* terms of (11) are then the downdated smoothing coefficients 
and replace the old smoothing coefficients. This process is repeated for j = k + 2, k + 6 , . . . , JS . ine 
terminal data equation yN = H n xn + v N is then downdated from the terminal SRIF array (3J. 


Implementation 


It is efficient to perform the process noise updates one component at at time. For example, in 
a four state filter with three process noise parameters, the smoothing coefficients from time update 
j + 1 can be expressed as 
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where the ()** terms are a reminder that all data has been processed in the SRIF. Note that the 
update order of the process noise parameters is pi, p2, p3. 


(5). 


The first step is to express the data equation to be removed in terms of using equation 
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The smoothing coefficients associated with pi are then downdated with an orthogonal transforma- 
tion. 
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For the sake of notation, the following variable (re)assignments are necessary: 
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Next, the smoothing coefficients associated with p2 are then downdated. 
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Finally, the smoothing coefficients associated with p3 are downdated. 
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Note that all smoothing coefficients before time update k + 1 remain unchanged. The ()* terms 
are the downdated smoothing coefficients and replace the old smoothing coefficients. The ensuing 
data equation y t+1 = H k+X z x+ , + I't+i is then used to continue this process to downdate the 
smoothing coefficients from time updates ib + 2 through N. In the end, the terminal data equation 
y N = H n x n + v N is then downdated from the terminal SRIF array [3]. 


Constructing the Orthogonal Transformation 


The orthogonal transformation in equations (17), (19), and (21) operates as 
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where the ()(Jb) terms are scalars and the ()* terms represent row vectors. As usual, let the 
orthogonally packed measurement be stored in the last column of the row vectors. To solve for the 
transformation T , multiply equation (22) by T ^ and take transposes. 
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Let A* = (A, a), where a = >/l - A 2 and A is a scalar to be determined. Note that the norm of A 
is 1. The transformation T is then constructed as an elementary Givens reflection. 

(24) 



With T constructed as such, equation (23) is then post-multiplied by A. 
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The scalar A and vector H can now be determined. 
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If instead of downdating, it is desirable to update the filter solution with an additional mea- 
surement, an elementary Givens reflection would be applied to the right side of equation (22). 
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Efficiency 


The discussed method was implemented on a VAX 8530. The filter test case required 30 minutes 
of CPU time to process 5800 data points in a 5.5 hour data arc of a typical GPS data processing 
scenario for the TOPEX/POSEIDON demonstration [6]; the filter estimated nearly 300 parameters. 
Downdating points from the resulting filter output required on average 10 seconds of CPU time per 
point. It proves that the linear combination of the data equation and the smoothing coefficients 
is relatively fast compared with downdating the terminal data equation from the terminal SRIF 
array. 

Numerical Integrity 

Stewart [5] has shown that downdating is stable in the presence of rounding errors. However, as 
he has also shown, if the spread of singular values is greater than half the computational precision, 
the precision of the smaller singular values may be lost. Fortunately, the downdating algorithm 
provides a way of detecting such ill-conditioned problems; the value of a (above and a similar 
quantity in the terminal downdating algorithm) is a reliable indication of trouble. Experience 
has shown that if a is less than 10 s , on a machine with a 15 digits of precision, the results 
of particular estimates may be inaccurate. This situation occurs when all information is removed 
from an estimated parameter. However, in an operational environment there are generally sufficient 
measurements or a priori information to avoid such situations. In these cases, experience has shown 
that the estimates are accurate to better than 10 digits. 

Applications 


In spacecraft orbit determination, often critical measurements, such as Very Long Baseline 
Interferometry (VLBI) measurements, are not available until long after the usual radio-metric data 
(e.g. doppler, range ) has been obtained. Using the discussed method, these critical measurements 
may be efficiently added after the usual data has been processed. When using optical data for orbit 
determination, camera pointing is often modeled as a white noise stochastic process. The discussed 
method permits the analyst to efficiently remove optical frames and replace them with others as 
desired. For large state systems with many stochastic parameters, such as GPS applications, outliers 

may be removed more efficiently and just as effectively by downdating rather than reprocessing the 
entire data set. 

Conclusion 

A new method to downdate a time- varying SRIF filter is presented. The method combines 
the algorithm of Gill [3] to downdate a matrix factorization and the decentralized SRIF algorithm 
of Bierman [4] to combine the results of independent time-varying SRIF filters. This method 
permits efficient removal of either measurements or a priori information. Additionally, updating 
a time-varying SRIF filter with either additional measurements or a priori information proceeds 
analogously. In both cases, a data equation is propagated by alternating between time mapping 
the equation and forming a particular linear combination with the SRIF’s smoothing coefficients. 
The terminal data equation that results can then be merged with or removed from the terminal 
SRIF array. Smoothing then proceeds as usual [2], For large state systems with many stochastic 
parameters, the discussed method is expected to be a critical component in the real-time reduction 
and analysis of large volumes of tracking data. 
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